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I Abstract 

■ 

We have developed spectral line On-The-Fly (OTF) observing mode for the Nobeyama Radio 
Observatory 45-m and the Atacama Submillimeter Telescope Experiment 10-m telescopes. Sets of digital 

1 autocorrelation spectrometers are available for OTF with heterodyne receivers mounted on the telescopes, 
including the focal-plane 5x5 array receiver, BEARS, on the 45-m. During OTF observations, the an- 

t-H ■ tenna is continuously driven to cover the mapped region rapidly, resulting in high observing efficiency and 

accuracy. Pointing of the antenna and readouts from the spectrometer are recorded as fast as 0.1 second. 
In this paper we report improvements made on software and instruments, requirements and optimization 
' of observing parameters, data reduction process, and verification of the system. It is confirmed that, using 

optimal parameters, the OTF is about twice as efficient as conventional position-switch observing method. 
• ■ Key words: radio lines: ISM — techniques: image processing — telescopes 

(N ■ 
t— i ; 

C^ - . 1. Introduction such as fast readout from the spectrometer and dealing 

with huge data, had prevented OTF from being applied 

J> , On-The-Fly (OTF) is a technique to perform mapping to spectral line observations. In the last decade spectral 

' observations efficiently with single-dish radio telescopes, line OTF has been made practicable in some observato- 

rS , In OTF, the antenna is driven continuously in a region ries (e.g., Mangum, Emerson, & Greisen 2000; Muders 

■ to be mapped and the data are taken in a short inter- et al. 2000; Ungerechts et al. 2000) and is established as 
val, instead of integrating at discrete positions on the sky an efficient method. Likewise, with the Nobeyama Radio 
(such a conventional "step-and- integrate" method is here- Observatory (NRO) 45-m and Atacama Submillimeter 
after referred to as position- switch [PSW] observation). Telescope Experiment (ASTE 1 : Ezawa et al. 2004; Kohno 
OTF observations have advantages over PSW modes as 2005) 10-m telescopes, continuum OTF observations have 
follows: (1) Observing efficiency improves since dead time already been available. We enabled the telescopes to be 
of the telescope is reduced and a large number 10 2 ) operated in OTF observing mode, as reported in this 
of on-source integrations are made per an emission-free paper. The OTF observing system of the telescopes 
reference integration; (2) The system variation (e.g., at- now works with an array of digital autocorrelation spec- 
mospheric conditions, receiver gain, pointing of the tele- trometers (MAC; Sorai et al. 2000). Heterodyne re- 
scope) less affects the map since the entire map can be ceivers equipped with the telescopes, including the 25- 
covered in a short period; and (3) Since the data are ac- beam focal plane array receiver (BEARS: Sunada et al. 
quired more frequently than the Nyquist sampling rate, 2000; Yamaguchi et al. 2000) on the 45-m telescope, can 
the spatial information is not lost. Comprehensive re- 

i j. • nrprj , l • l i j The ASTE project is driven by NRO, a branch of National 

view and discussion on ULb technique nave been made . ■ , A, r 

/r>rln>7\ Astronomical Observatory or Japan, in collaboration with 

by Mangum, Emerson, & Greisen (2007) . University of Chile, and Japanese institutes including University 

OTF technique has been commonly and widely used of Tokyo, Nagoya University, Osaka Prefecture University, 

in radio continuum observations. Technical limitations, Ibaraki University, and Kobe University. 
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be connected to the MAC and be available for OTF ob- 
servations. 

The capability of OTF to obtain a high-quality map ef- 
ficiently is important in various scientific fields. For exam- 
ple, investigating how the initial mass function (IMF) of 
stars is determined is one of the most important issues in 
astronomy. Recent studies have suggested that the stellar 
IMF is well related to the mass function (MF) of dense 
molecular cloud cores (e.g., Ikeda, Sunada, & Kitamura 
2007), which are sites of star formation. Therefore de- 
tecting a number of cores in various star forming regions 
is needed in order to study the relation between star form- 
ing activity and MF or physical/dynamical conditions of 
the cores. In this case sensitive unbiased mapping ob- 
servations of wide fields are required. In addition, OTF 
data are ideal to be combined with interferometric data, 
since they preserve spatial information in the (spatially) 
low frequency regime, which is lost via interferometric ob- 
servations. Thus, the fidelity of high resolution images of, 
e.g., molecular cloud cores or external galaxies, taken with 
interferometers are significantly improved by being com- 
bined with OTF data. The fidelity (preserving total flux) 
is essential to various kinds of studies. For instance, in or- 
der to investigate the evolution of the interstellar medium 
(formation and destruction of giant molecular cloud asso- 
ciations [GMAs]) across spiral arms in external galaxies, 
it is required to detect not only GMAs, which are spa- 
tially confined to spiral arms, but also diffuse emission in 
intcrarm regions. 

In order to make spectral line OTF observations avail- 
able with the NRO 45-m and the ASTE 10-m tele- 
scopes, various improvements have been applied to the 
control software system, COSMOS-3 (Morita ct al. 2003; 
Kamazaki et al. 2005), and related instruments. In par- 
ticular, fast and synchronized control of instruments and 
handling of huge amount of data are the difficulties, as 
discussed later. The way to determine optimal observing 
parameters is also complicated. 

In this paper we report practical information and im- 
plementation of the OTF observing mode for the tele- 
scopes. Basic concepts and parameters of OTF observa- 
tions are introduced in §2, followed by a brief description 
on controls of related instruments in §3. Requirements on 
observing parameters, procedures to estimate sensitivity, 
and how to maximize observing efficiency are explained in 
§4. We describe the data reduction process in §5. Desired 
characteristics of convolution functions to regrid the data 
onto a map are discussed and some appropriate functions 
are shown. In §6 we describe the measurements of driving 
performance of the antennas, comparison with PSW map, 
and verification of frequency-switch observing mode. 

2. Basic Observing Process and Parameters 

An OTF scan pattern is schematically illustrated in 
Fig. 1. A constant-speed raster scan with a single-beam 
receiver is assumed. Like usual PSW observations, the 
standard "chopper-wheel" technique (Penzias & Burrus 
1973; Ulich & Haas 1976) is employed to calibrate the 




Fig. 1. An example of scan pattern of OTF observations. 
The gray rectangle is a region to be mapped with a set of 
horizontal scans. The case of iV s S ca? = 2 (an OFF is inserted 
before every 2 ON scans) is shown. 

antenna temperature T£ in which the atmospheric and 
antenna losses are corrected. For the chopper-wheel cal- 
ibration, a hot load (R) at the ambient temperature and 
blank sky (SKY) are observed during the observation in 
an appropriate interval. Data at an emission-free refer- 
ence position (OFF) are taken before every on-source scan 
(ON, hereafter simply "scan") or set of several scans. The 
OFF is usually taken by pointing the antenna outside the 
source. However it can also be taken by shifting observing 
frequency (see §3.4.2). 

During the scan the antenna is driven at a constant 
speed (vscan) on the sky and the data are "dumped" from 
the spectrometer at an interval of idump- The "approach- 
run" at a speed v SC an is inserted before every scan to let the 
antenna move stably during the scan. If there is no OFF 
between two scans, a "transit-run" (from the end point of 
the scan to the start point of the next approach-run) is 
inserted. 

We set fundamental parameters as follows (see Fig. 1). 
In §4 we discuss how to determine them. The dimension 
of the mapping area is l\ x 1% (l\ along the scan, I2 across 
the scan). The scan speed of the antenna on the sky is 
Wscan- Time to be taken in an approach-run, a main (on- 
source) scan, and a transit-run are i app , t sca , a = ^i/v S can, 
and ttranj respectively. A separation between the scan 
rows is Al, the number of scans taken between a pair of 
OFFs is AfEQ, and the grid spacing of a map to be made 
is d. Time to be taken to slew the antenna between the 
mapped region and the OFF position is The pa- 

rameters i app , itran, and £°an depend on the performance 
of the antenna, the scan speed u SC an, position on the sky, 
etc., and arc typically in the range of a few to 10 seconds. 
Measurements of the performance of the antenna are de- 
scribed in §6.1. The time interval between each dump 
from the spectrometer is idump- 

3. Control of Related Instruments 

In OTF, spectral data are dumped at short interval (~ 
0.1 second) while the antenna is continuously driven across 
the sky. Consequently the following difficulties arise. 
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Fig. 2. Schematic illustration of three-layered COSMOS-3 
architecture. Rounded rectangles are components of 
COSMOS-3, cylinder symbols are data storage and files 
therein, and rectangles are individual instruments. 

Antenna. The antenna should be smoothly driven on a 
scan path. The sky coordinates, toward which the 
antenna pointed, have to be recorded for each data 
dump. 

Spectrometer. Data dump is done at short interval, and 
it must be synchronized with the antenna driving. 

Local oscillator. Dopplcr correction of the telescope 
with respect to a rest frame and frequency switch 
(FSW) arc difficult to be implemented in a conven- 
tional way, since exceedingly rapid control of the 
frequency is required. 

Data production and storage. High data production 
rate, caused by the fact that the data are dumped 
~ 10 2 times more frequently than typical PSW ob- 
servations, may cause trouble on data storage and 
reduction process. 

In this section we describe the control system and in- 
struments to resolve these issues. 

3.1. COSMOS-3 

The control software system of the telescopes is 
COSMOS-3 (Morita et al. 2003). It is designed to have 
a three-layer architecture, as illustrated in Fig. 2. The 
interface with observers is provided in the top layer. In 
the middle layer, MANAGER and MERGER exist. The 
MANAGER controls the observation by receiving instruc- 
tion from observers and interacting with Local Controllers 
(LCs), which directly control instruments, in the bottom 
layer. The MERGER receives outputs from LCs and pro- 
duces calibrated data. The observer can view a Quick 
Look (QLOOK) of the data and check if the observation 
is running properly. Synchronization of instruments with 
each other is managed with 1, 10, 50, and 100-Hz tim- 
ing pulses, which are generated and distributed from the 
standard clock system of the observatory. 



3.2. Antenna 

Before every approach-run, the antenna LC instructs 
the antenna to settle on the starting point of the approach- 
run. When the antenna is judged to be steady (the posi- 
tional error has been within a tolerance for a specified 
duration), the MANAGER decides the timing and co- 
ordinates of the following approach-run(s), scan(s), and 
transit-run(s). The antenna LC instructs the antenna to 
follow an ideal scan path at the decided timing. The ideal 
coordinates, toward which the antenna is desired to point, 
is hereafter called PROG values. The PROG values for 
every 0.1 second are calculated and transmitted to the 
antenna, 0.2 second before it is used. The antenna man- 
ages the timing to use the transmitted coordinates with 
the 10-Hz (0.1-second) timing signal. Two neighboring 
PROG values are linearly interpolated. The antenna is 
driven to follow the interpolated positions, referring the 
50 and 100-Hz timing signals. 

Meanwhile, the actual pointing of the antenna (here- 
after referred to as REAL values) jitters around the PROG 
values. The amount of the jitter (PROG — REAL) is mea- 
sured, in case of the 45-m telescope, using the encoder 
readout of the master collimator and the pointing devia- 
tion between the master collimator and the main reflector. 
In case of the 10-m telescope, the encoder readout of the 
antenna is used. The PROG (in both equatorial and hor- 
izontal coordinates) and PROG — REAL (horizontal coor- 
dinates) values are both written onto a file, the Antenna 
Log file, along with time stamps at every 0.1 second. 

3. 3. Spectrometer 

It is essential to synchronize the timing of data acquisi- 
tion with that of the antenna driving, particularly in case 
of OTF, since the data dump time is very short. The 
timing is managed as follows. 

When an observation starts, observing parameters, such 
as interval of data dumps (0.1 second) and number of 
datasets to be dumped, are calculated by the MANAGER 
from instructions made by the observer. They are sent to 
the spectrometer LC and then to the spectrometer, in ad- 
vance of data acquisition. For each scan, the MANAGER 
decides the time to start the scan as mentioned above. 
The LC is informed of the time, waits until 1 second be- 
fore it, and then sends a command to the spectrometer 
to start integration. After receiving it, the spectrometer 
begins to acquire the data, triggered by the next fall-down 
of the 1-Hz, 50%-duty timing pulse. It dumps the data at 
every 0.1 second during the scan. 

Though the spectrometer outputs datasets at every 0.1 
second, observers can specify idump to be O.liV (N = 
1,2,...) seconds to reduce the data size. Every 0.1-second 
output is transferred to the LC. The LC averages succes- 
sive N datasets if N > 1, and writes them down onto a 
file (the Spectrometer Log file) along with time stamps. 
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3-4- Local Oscillator 

3.4-1- Doppler Tracking 

When a wide area is mapped, the radial velocity of the 
telescope (v iaJ d) with respect to the reference frame (e.g., 
the Local Standard of Rest [LSR]) significantly changes 
from point to point: the velocity gradient across the sky 
plane amounts to roughly 0.8 km s _1 deg -1 (LSR) at max- 
imum in NRO. In PSW observations, the radial velocity of 
the telescope at each position is calculated and the local 
oscillator (LO) frequency is shifted before the integration 
starts, in order to track the Doppler shift caused by Av ta/ &. 

However, in OTF observations, the antenna runs on a 
scan path (thus t> ra d gradually changes), and during which 
the spectral data are continuously taken at a short inter- 
val. In this case it is difficult to track v ra a by shifting the 
LO frequency continuously during the scan. Therefore 
we do not make any shift on LO frequency during the 
OTF scans. Instead, the Doppler correction is carried out 
by software: in the production process of spectral data 
(MERGER; see below), w rac j for each ON integration is 
calculated, and then channel shift operation is made. 
3-4-2. Frequency Switching 

When spatially-widespread and/or narrow lines (e.g., 
CO in nearby molecular clouds) are mapped, FSW ob- 
serving technique has been widely used. Conventional 
FSW observations with the NRO 45-m and ASTE 10-m 
telescopes are performed as follows. The LO frequency 
is switched between the original value and slightly 
shifted value v^ ± Av-^o . The switching cycle is in the 
order of ~ 1 Hz, and the frequency throw Av^o is typically 
several MHz. The averager/integrator of the spectrom- 
eter behaves in adding and subtracting mode, synchro- 
nized with the frequency shift of the LO, during the ON 
integration (the averager/integrator does nothing during 
a margin for the transition of frequency). Thus the ON 
concurrently plays a role of OFF. Accordingly the origi- 
nal signal and a negative-amplitude signal appear on the 
spectrum at a separation of ZWlo- The final spectrum is 
demodulated by invert, shift and add operation. 

In OTF, the difficulty of shifting LO frequency during 
a short integration again arises. We have implemented 
FSW observing mode in OTF by shifting the LO fre- 
quency at OFF position only: v^ is used during ON, 
while v®q ± Av\^o is used during OFF. The emission line 
existing (if any) at OFF position shifts by Av^o, thus it 
can be virtually regarded as emission-free reference. 

3.5. Spectral Data Production 

As described above, the antenna pointings are recorded 
in the Antenna Log file at every 0.1 second, and the 
readouts from the spectrometer are written onto the 
Spectrometer Log file at every O.liV seconds. The 
MERGER, which is a process to produce spectral data 
files, runs during the observation (it can also be executed 
asynchronously after the observation). The MERGER 
process reads the Spectrometer Log file, calculates the 
standard chopper-wheel equation to calibrate the antenna 
temperature, and shifts the resultant spectra along the fre- 



quency channels to correct Doppler shift. The MERGER 
also reads the Antenna Log file. By comparing time 
stamps written in both log files, it calculates sky coor- 
dinates, toward which the antenna was actually pointed 
at the moment of each dump, using the PROG (equato- 
rial) and PROG — REAL (horizontal) values. The derived 
coordinates are written in the header of the spectral data 
file. In this way intensity- and velocity-calibrated spectra 
are written on a file, called the Raw Data file. 

In practical implementation, the trouble with the Raw 
Data is their large sizes. If the 25-beam, 1024-channel 
data are represented as 32-bit and dumped every 0.1 sec- 
ond, the data rate goes up to ~ 3.5 GBhour -1 (data 
headers are neglected). In order to reduce the size of 
the Raw Data, the quantization level, Q [bit], must be 
as small as possible. When the peak antenna temper- 
ature of the intrinsic spectrum is T p and rms noise is 
a, observed antenna temperature between the minimum 
(— TVcr; N ~ 3) and maximum (T p + Na) is quantized into 
2*3 levels. The quantization noise is roughly written as 
o~q = (T p + 2Na)/2^ +1 . Since the noise level increases to 

yV 2 + Oq. the overhead in the observing time (the ratio 

between Q-bit and oo-bit cases) is 

tq_ . (T n 



1 



-r2Na\ 2 
2Q+ 1 a J ' 



(1) 



meaning that the loss increases when the signal-to-noisc 
ratio of individual spectrum, T p /a, is high. Supposing an 
extreme case for the current system: peak antenna tem- 
perature T p = 100 [K], system noise temperature T sys = 
100 [K], frequency resolution B=l [MHz], and data dump 
time idump = 0.1 [s], we obtain ig/ioo - 1 ^ (320/2 Q+1 ) 2 . 
For Q = 8, 12, and 16, tq/t^ — 1 becomes 0.40, 1.5 x 
10 ~ 3 , and 6.1 x 10~ 6 , respectively. We adopted Q = 12 
(4096 levels), for which the quantization noise is practi- 
cally negligible, resulting in the maximum data rate of 
~ 1.3 GB hour Jl . 

In addition, the MERGER supports options of channel 
trimming and channel binning. When these options are 
specified by the observer, the MERGER trims and/or bins 
up channels of the spectra before writing them down on 
the Raw Data file. 

It should be noted that the 12-bit quantization may 
affect the data in the following situation: bandpass be- 
comes nearly (i.e., diverges) at the band edges, or an 
extremely intense spurious signal appears. The observer 
should pay attention to the obtained spectra through 
QLOOK. This problem can be evaded by trimming the 
band edges and/or spurious signals, since the data are 
quantized after the trim. 

4. Observing Parameters 

4-1- Requirements on Sampling 

As a result of an observation, the mapped region is 
filled with data points. The data sampling separation is 
fscan^dump along the scan, Al across the scan. Mangum, 
Emerson, & Greisen (2007) discussed requirements on 
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the sampling, as summarized as follows. It is requested 
that, at least, the sampling rates fscanidump and Al are 
both more frequent than the Nyquist sampling rate X/2D, 
where A is the observed wavelength and D is the diameter 
of the antenna aperture. In case of A = 2.6 mm (115 GHz) 
observations with the 45-m telescope X/2D ~ 6", which 
corresponds to ~ 1/2.5 of the half-power beam width 
(HPBW). Practically a sampling more frequent than the 
Nyquist rate is required to avoid aliasing noise and beam 
smearing effects. Since the data points do not align on 
any regular grid due to the antenna jitter etc., the data 
should be rcgridded onto a regular grid using a gridding 
convolution function (GCF) in the data reduction process 
(see §5). 

4-2. Estimation of Sensitivity 

We estimate the sensitivity of a single-beam observa- 
tion. Application to a multi-beam receiver is discussed 
later. 

The total number of scan rows in an observation is 

A rc 

becomes 



1%/Al + l. The total on-source integration time 



ON 



AW 



ow^scan- 



(2) 

The total time spent to run an observation including R, 
SKY, OFF, antenna slew, etc. is estimated to be 



The factor 77 is a constant determined by the extent of 
the used GCF and is calculated as follows. Suppose 
that observed points i — 1,2,... are uniformly distributed 
around the grid point and each point has a spectrum 
Ti(k) (k = I, Nch), rms noise temperature o~i, and a 
GCF weight m;. We assume that the on-source integra- 
tion time to and therefore the rms noise temperature er.; = 
do = Tgy S /(r]q-\/Bto) of each point are both constant. Here 
T sys is the system noise temperature, B is the frequency 
resolution of the spectra, and ry q is the quantization effi- 
ciency of the spectrometer. In case of MAC r/ q — 0.88; and 
hereafter i] q is omitted from expressions. The convolved 
spectrum T(k) is written as T = Q2 w iTi)/Q2wi), and 
its noise temperature a becomes a = (\/J2 w i /J2 w i) a o = 

WV^SK' whcrc *SS = *o (E^)7£K 2 )- If we take 
the grid spacing as the unit of spatial length and re- 
define <o a s the on-source integration time per unit area 
(1 grid cell), summations can be rewritten with integrals: 
^ccU = ^0 (/ wdxdy) 2 1 J w 2 dxdy = rjto- Approximate values 
of r] for GCFs Besselx Gauss, Sincx Gauss, Gauss, Pillbox, 
and Spheroidal (see §5) with default parameters are, re- 
spectively, 4.3, 1.2, 6.3, 1.0, and 10. 

Redefining B as the frequency resolution of a map to be 
made, the noise of the map due to on-source integration 
is estimated to be 



OBS 



■ ^OH 



N, 



OFF 
SEQ 



/c; 



(3) 



ATI (0 N ) 



sys 



(9) 



where toFF is an integration time for an OFF, and / ca i 
is an overhead of R-SKY calibration (if 1 minute is con- 
sumed to obtain R and SKY data at every 15 minutes, 
/cai = 16/15). The ton is an overhead time per one scan 
row, which consists of go-and-return to the OFF point 
2^tran j time for approach- and transit-run i app and ttran, 
thus is written as 



the standard radiometer equation. On the other hand, the 
number of OFF points used to consist a map grid point 
is roughly written as dj Al (here the extent of the GCF is 
neglected). Thus the effective OFF integration time for a 
grid cell is 



^oh — 



2t 



OFF 
trail 



N. 



SEQ 



ATSEQ _ x 



,OFF 
r cell 



Al 



*OFF, 



app - 



Af 



SEQ 



-t. 



(4) and the noise due to OFF points becomes 



The t^an depends on the distance between the mapped Z\T£(OFF) = 
region and the OFF position; i app and t tran should be 



chosen so that the antenna follows the PROG position 
during the scans. Now the ratio of on-source time to the 
total time spent is 



(10) 



(11) 



Therefore the total noise level of the map is written as 



770N/OBS 



t 



OX 



ATI = \M7a( on ) 2 + ^Ta(OFF) 2 



fOBS 



tsc 



rSEQ 

v scan 



1 

fcal 



(5) 
(6) 



T / 1 1 

sys 



(12) 
(13) 



£oh + £off/A^ 

The total on-source integration time for a map grid 
point, tJSgjj, is a sum of time during which the beam scans 
within the grid cell. Since the data are convolved using a 
GCF to construct a regularly gridded map, effectively a 
factor 77 is multiplied: 

d 2 



fu\ t ON +OFF ■ 
V-D y coll ''coll 

Validity of this equation has been practically confirmed 
(see §6.2). 

4-3. Optimization 

The noise level of a map achieved in unit observing time 
AT^(0) is written as 



ON 
cell 



-t 



ON 



(7) 



ATl(0)=ATiJt°™ 



(14) 



.sys 



1 1 
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Fig. 3. Dependence of AT^(0) (normalized by the mini- 
mum) on toFF are shown in three CcLSCS: Cscan = 20 (solid), 40 
(dashed), and 60 (dotted) seconds. Common parameters are: 
Z 1= 600 ["], Al = 5["], d = 7.5["], t OH = 25 [s], and N^ = l. 



Bcsselx Gauss convolution (r) = 4.3) is assumed. 



i it i toFF 
! -scan i 6 OH i SPO I iv iowJcal 



N SEQ 

y \ jv scan 

and is minimized when ^off is optimal: 
-?—ATl(0) = 

OTOFF 

leads to 



t: 



optimal 
OFF : 



(^scan + *Oh) 



7}dt s 



h 



SEQ 



(15) 



(16) 



(17) 



This formula is a generalization of the well-known relation 
ioFF = VA^on for PSW observations, where N is the 
number of ONs taken per one OFF. Toward toH — * 0, Eq. 
(17) resolves itself into the viV-relation. 

Dependence of AT^(0) on toFF is plotted in Fig. 3. 
For cases of i scan = 20, 40, and 60 seconds (other param- 
eters are shown in the caption), Z\T£(0) is minimized at 
7, 12, and 17 seconds, respectively. If ioFF is 



optimal 
OFF 



t 

shorter than the optimal value, ZiTJ(OFF) dominates the 
map. On the other hand, <off longer than toF™^ ^ s 
excessive since the noise level of the map is limited by 
AT£(ON). 

Using the above notations, it is quantitatively shown 
how the observing efficiency improves, by adopting ap- 
propriate parameters, compared with PSW observations. 
The following two factors contribute: (1) The ratio of 
on-source time to the total time spent, t?on/obS: t> e ~ 
comes larger because not only the dead time (antenna 
slew, etc.) is reduced, but also the OFF integration time 
is relatively shorter (£off -C t scan ); and (2) In general 
*?en < *?eU F ; thus AT A ncarl y e q uals to ^a(ON) (in- 
stead of ATI - \/2/ATX(ON), which is applicable to PSW 
observation with ton = ^off) ■ As ton — ► and < SC an — * oo, 
the both factors respectively correspond to reduction of 
observing time by a factor of 2. Thus OTF is, theoreti- 
cally, up to 4 times more efficient than PSW. In practice, 
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Fig. 4. A schematic illustration of OTF scan patterns with 
BEARS. The array is inclined by an angle 8 with respect to 
the direction of scans in order to make the scans by neighbor- 
ing beams separated by Al. The first scan (labeled as "Scan 
1") and the next ("Scan 2") are offset by 5AI (see text). 

improvement of efficiency amounts to a factor of ~ 2. 
4-4- Application to an Array Receiver 

The above discussion is made for a single-beam receiver. 
In case of an array (multi-beam) receiver, some expres- 
sions change. Here we consider the case of BEARS, a 
5x5 focal plane array. 

Fig. 4 schematically shows OTF scans with BEARS 
considered here. The array is inclined with respect to 
the scan direction by an angle 9: a neighboring couple 
of beams makes a pair of scans separated by a distance 
Al = Ls'mO (L = 41"1 is a beam separation). The next 
scan runs 5AI away, thus the number of scan rows is writ- 
ten as A row = l2/(5Al) + 1. In this case scans made by 5 
beams in each row of the array fill the mapped region at a 
separation of Al. Since the 5 rows of the array respectively 
cover the map, 



ON 



hAl 



and 



t 



OFF 



5d 
Al 



ioFF- 



(19) 



The noise level of the map is obtained by substituting tx-jj 
and tg FF in Eq. (13) with Eqs. (18) and (19). Following 
the transformation in the previous subsection, we have 



t 



optimal 
OFF 



(£scan +^Oh) 



7]dt 



I tvtSEQ 



(20) 



the same expression as Eq. (17). 

The strategy mentioned above focuses on mapping 
speed: compared with the case of a single-beam receiver, 
the mapped region can be covered 5 times faster and the 
integration becomes 5 times deeper. There may be an- 
other strategy, conscious of the uniformity of the map. 
Using the same parameters as the case of a single-beam 
receiver (each beam runs at a separation of Al), the in- 
tegration simply becomes 25 times deeper in the same 
observing time. The point is that every scan path is ap- 
proximately traced by all the 25 beams: characteristics 
of the beams are expected to be averaged out. However 
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whether this advantage is realized or not depends on the 
condition, since the system variation may spoil the uni- 
formity if the observing time becomes too long. 

5. Data Processing 

5.1. Data Reduction 

Reduction process of OTF data is done with NOSTAR 
(Nobeyama OTF Software Tools for Analysis and 
Reduction) developed in NRO. It is designed to run on 
UNIX or UNIX-like operating systems. Its core func- 
tions (baseline subtraction, creating cube FITS from spec- 
tra, etc.) are provided as command- line tools written in 
C /FORTRAN in order to enable batch processing. These 
command-line tools are wrapped in graphical user in- 
terface (GUI) written in the Interactive Data Language 
(IDL). 

A run of observation produces a Raw Data file (see 
§3). Spectra from all the used spectrometers are contained 
within it. The first step of data reduction process is to ex- 
tract the data to the user's working directory. The task, 
named Split, does not simply copy the Raw Data file, but 
divide it into Split Raw Data files according to the spec- 
trometers. Namely, a Raw Data taken with BEARS is 
split into 25 Split Raw Data files. Accordingly the follow- 
ing process can be separately applied for each Split Raw 
Data. 

Subsequent procedures of data reduction (baseline sub- 
traction, bad data flagging, etc.) are basically the same as 
those for PSW data, except for large size of the data. The 
GUI is designed to be batch-oriented, to help users to pro- 
ceed the reduction process rapidly without being bothered 
with a large quantity of data. Each process overwrites the 
Split Raw Data itself in order to avoid running out of space 
of the working directory 

Finally a map (cube FITS) is made from the processed 
spectra. An OTF observation (or a series of observations) 
produces a set of data points which fills the mapped region 
with spacings smaller than the Nyquist sampling rate, as 
described in §4. The data are convolved into the map 
using a GCF. Desired characteristics of GCFs and appro- 
priate function forms are discussed below. The obtained 
map may suffer from the so-called scanning noise along 
the scan direction, in addition to the statistical noise. The 
scanning noise can be effectively removed by combining 
two maps made from orthogonal scans using the so-called 
basket-weave method. We have implemented the PLAIT 
algorithm described by Emerson & Grave (1988). 

5.2. Gridding Convolution Functions 

GCFs used to make maps are desired to have the fol- 
lowing characteristics. 

First: form of GCF is similar to that of the telescope 
beam itself. Convolution with such a GCF corresponds to, 
in the Fourier domain, that a weighting function (Fourier- 
transformed beam pattern) is multiplied twice to the in- 
trinsic spatial frequency distribution. Consequently the 
best signal-to-noise ratio is achieved. Though we cannot 
know the beam pattern at infinite accuracy in practice, 



one should choose a GCF so that it mimics the telescope 
beam. 

Second: the GCF's energy concentration ratio is high. 
In practical convolution operation, the extent of the GCF 
is finite. Thus, in the Fourier domain, the GCF has artifi- 
cial frequency components (sidelobes). By regridding the 
data, these sidelobes are folded onto the primary compo- 
nent, resulting in the so-called aliasing noise. In order to 
preserve the observed spatial frequency information, the 
aliasing effect must be as small as possible. As an index, 
energy concentration ratio 

f A \C(r))\ 2 dv 



irjc( V wd V 



(21) 



is introduced (Briggs, Sramek, & Schwab 1999), where the 
GCF is c(l), I = x/Ax, Ax is the grid spacing, C(rj) is the 
Fourier transform of c(l), and A is the area in which r\ < 1. 
Here TZ represents the degree of concentration of the GCF 
within A: it is expected that the larger TZ is, the smaller 
the aliasing effect is. 

The following GCFs (their shapes are shown in Fig. 5) 
arc implemented so far: 

Besselx Gauss: a Gaussian-tapered Jinc function 



c(r): 



Jjjjrr/a) 
Trr I a 





exp 



-(£) 2 J (r<iU) 
(otherwise) 



(22) 



where J\ is the lst-order Bessel function, r is the 
distance between the data point and the grid point 
(the unit is the grid spacing). Parameters a, &, and 
i? max can be arbitrary chosen, and are set as a = 
1.55, b = 2.52, i? max = 3 by default (see below). 
Sincx Gauss: a Gaussian-tapered Sine function 



c(r): 



sin(7rr/ a) 
7rr I a 

o 



exp 



(r < -Rmax) 

(otherwise) 



.(23) 



The default values for parameters a, 6, and i? ma x are 
the same as those for the Besselx Gauss function. 
Gauss: a pure Gaussian 



c(r) 




The default parameters are a 
Pillbox: a cell-averaging function 



(r < i?max) 

(otherwise) 

1 i ^ma: 



c(x,y) 



1 (M<0.5,|y|<0.5) 
(otherwise) 



(24) 



(25) 



where x and y are the distance between the data 
point and the grid point (the unit is the grid spacing) 
in the Cartesian coordinates along the map grid. 
Spheroidal: the spheroidal function 

|l-^| Q V>ao(c,7fe)x 
|l-^| a- 0ao(c,?7y) 

(|%|<i,kl<i) 

(otherwise) 

where rj x = 2x/m, rj y = 2y/m. The parameters m(= 
4,5,6,7,8) and a(= 0.0,0.5,1.0,1.5,2.0) define the 



c(x,y) 



(26) 
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Fig. 5. Shapes of implemented GCFs c(x): Besselx Gauss 
(solid), SincxGauss (dashed), Gauss (short- dashed), Pillbox 
(dotted), and Spheroidal (dot-dashed); normalized so that 
c(0) = 1. Function parameters are their default ones (see text). 

shape of the function. See Schwab (1984) for details. 
The defaults are m = 6 and a = 1.0. 

The first four GCFs (Besselx Gauss, SincxGauss, 
Gauss, and Pillbox) arc taken from Mangum, Emerson, 
& Greisen (2007). Given the cut-off radius i? max = 3, the 
default parameters a = 1.55, b = 2.52 for the SincxGauss 
function are derived so that the energy concentration ra- 
tio 1Z becomes the maximum (Schwab 1984). Using the 
same parameters, 1Z does not become the very maximum 
for the Besselx Gauss function (the maximum is found at 
a ~ 1.4, b ~ 2.7), but is still high enough. 

The spheroidal function is taken from Schwab (1984), 
and is often used to regrid the interferometric visibility (u- 
v) data. The function form is derived so that 1Z becomes 
as high as possible. 

As a result of convolution, spatial resolution of the map 
becomes lower than the telescope beam. We estimate the 
influence of the beam broadening by calculating the re- 
sponse to a point source. Fig. 6 shows the peak temper- 
ature and full-width at half-maximum (FWHM) of the 
source (i.e., effective beam) on the obtained map, as func- 
tions of the grid spacing d. If d is too small, effective 
integration time for a grid (t^u) becomes small, which 
leads to a large noise level. In this case the map is too 
much oversampled, since the spatial resolution is limited 
by the telescope beam. On the other hand, if d is too 
large, the effective beam broadens to ~ 2d. 



6. Verification 

6.1. Performance of the Antennas 
In 62, 



and tt 



we introduced two parameters t app aim Hran, 
which are duration of "approach-run" and "transit-run" , 
respectively. They depend on the driving speed and 
performance of the antenna. In this subsection we de- 
scribe the measurements to determine t app and t tIan . The 
amount of pointing jitter (PROG — REAL) is also mea- 
sured. 



Som 


0.8 






Ph 


0.6 











-i 


0.4 


1 






0.2 


Oh 







Peak Temperature 




Effective Beam 



















2.5 



2.0 | 



L.5 



1.0 



0.2 0.4 0.6 OS 1 

(grid spacing)/(beam FWHM) 

Fig. 6. Dependence of peak temperature of a point source 
(left-hand axis) and FWHM of effective beam (right-hand 
axis) on the grid spacing. Beam smearing effect is neglected. 
The vertical axes are normalized by the original (no convo- 
lution) values. The telescope beam is assumed to be a pure 
Gaussian, and the Besselx Gauss convolution with default pa- 
rameters a = 1.55, b = 2.52, and i? m ax = 3 is used. 
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Fig. 7. The PROG - REAL difference of the 45-m telescope 
along the scan direction on the sky is shown against the time 
elapsed after starting the approach-run. The measurements 
were done at every 0.1 second. 

6.1.1. NRO 45-m 

First, we made scans along both Azimuth (Az) and 
Elevation (El) at various scan speed v scan , in order to 
determine t app . The driving speed of the antenna is 
u S can/cos(El) for Az scans, and u scan for El scans. Fig. 7 
shows some examples. The PROG — REAL difference on 
the sky is plotted against the time elapsed after starting 
the approach-run for three cases (a) El scan with i> scan = 
40 ["/s] at El = 54°; (b) El scan with u scan = 160 ["/s] at 
El = 54°; and (c) Az scan with u scan = 240 ["/s] at El = 74° 
(the driving speed is 870"/s along Az). The case (a) rep- 
resents the typical scan at 115 GHz, while the case (c) 
corresponds to the maximum scan speed at 22 GHz. The 
antenna runs stably after the initial delay and (for large 
driving speed) an overshoot. The time spent until the sta- 
ble run begins is adopted to be t app . Determined t app for 
the three cases are, respectively, 5, 7, and 11 seconds. 
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Fig. 8. Derived t a pp and ttran for the 45-m telescope. The 
horizontal axis shows the driving speed, u SC an/cos(El) for Az 
scans and « scan for El scans. 

Secondly, itran is determined in a similar way as t app : 
back-and-forth scans along Az and El are made using ap- 
propriate t app . Various tt ran are tried. The optimal itran 
is defined as the shortest one for which the antenna starts 
to move stably after the next approach-run. Fig. 8 shows 
derived i app and itran- The both of them monotonically 
increase with the driving speed. A slight jump of f app is 
found at the driving speed of ~ 200 ["/s\. It is due to 
"slow-start slow-stop" control of the telescope, which is 
implemented not to give sudden and large acceleration to 
the antenna. 

It is found that the PROG — REAL value does not con- 
verge into when the driving speed is large (see Fig. 7). 
Instead, PROG — REAL offset becomes almost constant, 
which amounts to |PROG-REAL| ~ 0.016 [s] x v scan . The 
offset does not affect the observations, since it is about 1/6 
of the sampling separation along the scan for £dump of 0.1 
second. Excepting the offset, the jitter of the antenna 
pointing is within a few arcseconds. The error of the sam- 
pling separation is determined by the differential of the 
jitter between the neighboring sample, which is almost 
within l"-2". 
6.1.2. ASTE 10-m 

For the ASTE 10-m telescope, t app , itran, and the point- 
ing jitter are measured in a similar way as the 45-m. The 
iapp measurements for the cases (a) El scan with v scan = 
50 ["/a] at El = 30° and (b) Az scan with u scan = 100 ["/s] at 
El = 70° (the driving speed is 290"/s along Az) are shown 
in Fig. 9. The case (b) corresponds to the maximum scan 
speed at 350 GHz. The PROG - REAL values converge 
into ~ in a few seconds, with jitters of < ±1". From 
this result, we adopted 4 seconds as i app . The ttran is also 
measured: for practical observing parameters, 2 seconds 
are enough. If receivers for higher frequency is installed in 
the future, the antenna performance should be measured 
again. 



Fig. 9. Same as Fig. 7, but for the 10-m telescope. 

6.2. Comparison with PSW Map 

In order to confirm the validity and efficiency of OTF 
observations, the same field has been observed with 
both OTF and PSW, and the resultant maps are com- 
pared. Here we describe ASTE CO J = 3-2 (345.8 
GHz) test observations toward a field centered at (l,b) = 
(37 o 45',-0°12'). 

The observations were made in 2005 August and 
September, as a part of a CO J = 3 — 2 Galactic plane 
survey (Sawada et al., in preparation). The frontend was 
a cooled SIS mixer receiver, SC345. The system noise 
temperature was typically 200-300 K in double sideband 
(DSB) during the observations. The HPBW of the tele- 
scope at 350 GHz was measured to be 22". The main 
beam efficiency was « 0.6. The backend was a 1024- 
channel MAC, which covers an instantaneous bandwidth 
of 512 MHz (440 km s -1 ) with a spectral resolution of 
1.0 MHz (0.87kms _1 ). The pointing of the telescope 
was calibrated by tracking a compact CO source W Aql 
in every 1 or 2 hours, and was within the accuracy of 
5". We performed intensity calibration by observing a 
standard source, Ml 7 SW, in every 2 hours. The re- 
producibility of T£ was 5% (la). By comparing the 
observed spectra with those measured with the Caltech 
Submillimetcr Observatory (CSO) 10.4-m telescope with 
a single-sideband (SSB) filter (Wang et al. 1994), we ob- 
tained scaling factors to convert ASTE TJ(DSB) into CSO 
TX(SSB). Hereafter is shown in SSB scale. 

The OTF observations were carried out with parame- 
ters u 8C an = 50"/s and Al = 8". Two longitudinal scans 
and two latitudinal scans (in the Galactic coordinates) 
were made. We made a map, whose grid spacing is 
8" x 8" x 1 kms -1 , using the BesselxGauss convolution; 
and then the map was resampled onto a 10" x 10" grid to 
match the PSW data. A small portion of the region was 
observed with PSW: 9x9 points separated by 10". For 
each point 10-second integration was made twice. We have 
convolved the PSW data to make a map having the same 
resolution (25") and grid spacing as the OTF map. Since 
the convolution is highly incomplete at the outermost grid 
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Fig. 10. Velocity channel maps obtained from OTF (left) 
and PSW (right) observations. The velocity range of 
55.5 < «lsr < 61.5 km s _1 is sliced into Av = 1 km s — 1 chan- 
nels. Contours are Tj£ = 1 to 9 K at an interval of 1 K. 
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Fig. 11. A correlation plot between PSW and OTF T£ in 
the velocity range 50 < i>lsr ^ 70 kms . The error of 
±3<t is shown in bottom right. Solid line shows the best fit 
T*(OTF) = (1.045 ±0.004)T* (PSW). 

points, we use the inner 7x7 pixels for comparison. 

The la noise levels are 0.25 K (OTF) and 0.12 K 
(PSW). The obtained OTF noise level agrees with the 
one derived using the equations in §4, 0.26 K. It is proven 
that the system has achieved expected observing effi- 
ciency. Velocity channel maps of OTF and PSW are 
shown in Fig. 10. The OTF map is consistent with the 
PSW one. Fig. 11 shows a pixel-to-pixel correlation plot 
between them. A least-square fit gives the correlation 
TX(OTF) = (1.045 ±0.004)T£(PSW). It is confirmed that 
the OTF map agrees with the PSW map within the accu- 
racy of relative intensity calibration, 5% (reproducibility 
of the intensity of the standard source). 

6.3. FSW Observations 

We have implemented FSW observing mode in OTF as 
described in §3.4.2. As a test of FSW-OTF observations, 



IRAS 04369+2539 (IC 2087) in Taurus molecular cloud 
was observed. High-velocity wing (outflow) emission was 
found around the source in the 12 CO J = 1-0 line (Heyer et 
al. 1987). Mapping wing emission in widespread molecular 
cloud is a presumable science case for FSW observations. 
Stable spectral baseline is required. 

Observations were carried out with the 45-m tele- 
scope and BEARS. The HPBW is 15", and the sys- 
tem noise temperature was typically 350 K (DSB). The 
MAC was used in high- resolution mode, i.e., having 32 
MHz (83 kms -1 ) instantaneous bandwidth and 63 kHz 
(0.16 kms -1 ) resolution. The OFFs arc taken at start- 
ing point of approach-runs in order to reduce dead time 
to slew the antenna. The frequency throw of the LO, 
Avlo, was set to 12 MHz (corresponding to 31 kms -1 ). 
The TX(DSB) was converted into TX(SSB) by comparing 
spectra of a standard source measured with BEARS and 
an SSB receiver, S100. Linear baselines were subtracted, 
and a 7."5-grid map was made using the Besselx Gauss 
convolution. 

Fig. 12 shows a set of velocity channel maps. Redshifted 
wing emission toward north-south direction is success- 
fully detected beyond the ambient cloud velocity, wlsr — 
6 kms -1 . Fig. 13 shows line profiles at two positions 
(Aa,AS) = (-4', +6') and (-4',-T). At (-4', +6'), wing 
emission extending to ~ 14 kms -1 is seen, which is consis- 
tent with Heyer's results. On the other hand the profile at 
(— 4', — 1'), which sharply truncates at ~9kms , demon- 
strates a straight baseline in a velocity range of 40 km s -1 . 

6.4- Application to Weak/ Broad Line Objects 

In general, OTF is effective in particular when widely 
distributed and intense line is mapped, since a beam 
runs across a map grid within very short duration, typi- 
cally (a few-10) x 0.1 seconds, without any overhead to 
point discrete positions. However, application to rel- 
atively small-field (~ a few arcminutes square), weak- 
and/or broad-line (~ several 10 mK, several 100 kms -1 ) 
sources (e.g., external galaxies) has been successful. For 
example, Hirota et al. (2008, to be submitted in this vol- 
ume) observed a galaxy IC 342 using the NRO 45-m and 
BEARS, and achieved 8 mK rms in T£ at a velocity res- 
olution of 5 km s -1 . 

7. Summary 

We have made spectral line OTF observations available 
at NRO 45-m and ASTE 10-m telescopes. Digital auto- 
correlation spectrometers can be operated in OTF mode 
(the data sampling interval is as fast as 0.1 second) with 
heterodyne receivers mounted on the telescopes, including 
the 25-beam array receiver, BEARS. Improvements of the 
software and instruments to enable fast and synchronized 
controls (e.g., antenna driving, data acquisition, Doppler 
tracking, frequency switching) were described. Sensitivity 
of the obtained map was expressed using observing pa- 
rameters, and we showed how to determine and optimize 
the parameters. Performance of the antennas was mea- 
sured and was proven to be high enough for up to 115 
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GHz (45-m) or 350 GHz (10-m) observations. The OTF 
system has improved observing efficiency by a factor of 
~ 2 compared with PSW: its mapping capability opens a 
prospect in various fields of study. 

We thank NRO staffs and the ASTE team members 
for their support. In particular, we acknowledge Satoshi 
Hongo, Aya Higuchi, Rie Miura, Nobuyuki Yamaguchi, 
and Kunihiko Tanaka for their contribution to develop- 
ment and evaluation. 
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Fig. 13. The CO J = 1-0 line profiles at two posi- 
tions (Aa, AS) = (-4', +6') [top] and (-4', -1') [bot- 
tom, —2 K offset]. The positions are relative to 
(a,5)i950 = (4 h 36 m 54?6,25°39'17") 



